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Direct Numerical Simulations established that the FENE-P model of viscoelastic flows exhibits 
the phenomenon of turbulent drag reduction which is caused in experiments by dilute polymeric 
additives. To gain analytic understanding of the phenomenon we introduce in this Letter a simple 
1-dimensional model of the FENE-P equations. We demonstrate drag reduction in the simple model, 
and explain analytically the main observations which include (i) reduction of velocity gradients for 
fixed throughput and (ii) increase of throughput for fixed dissipation. 
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The addition of few tens of parts per million (by 
weight) of long-chain polymers to turbulent fluids can 
bring about a reduction of the friction drag by up to 
80% jjj. This "drag reduction" phenomenon has im- 
portant practical implications besides being interesting 
from the fundamental point of view, integrating turbu- 
lence research with polymer physics. In spite of intense 
interest for an extended period of time [0, ||, [| , Sreeni- 
vasan and White |Q recently concluded that "it is fair to 
say that the extensive - and continuing - activity has not 
produced a firm grasp of the mechanisms of drag reduc- 
tion" . Recently however it was shown that drag reduc- 
tion is observed in Direct Numerical Simulation of model 
viscoelastic hydrodynamic equations || ||, [jj . From the 
theoretical viewpoint these observations are crucial, indi- 
cating that the phenomenon is included in the solutions 
of the model equations. Understanding drag reduction 
then becomes a usual challenge of theoretical physics. 
In this Letter we present a further simplification of the 
model equations and gain analytic insights into the phe- 
nomenon. The FENE-P equation for the fluid velocity 
u(r, t) contains an additional stress tensor related to the 
polymer: 



du 
~dt 



+ (u-V)u = -Vp + v s V 2 u + V T + F , (1) 



where v s is the viscosity of the neat fluid, F is the forcing 
and the stress tensor T is determined by the polymer 
conformation tensor R according to 



f{r,t) 



Pi 



R(r,t) 



(2) 



here v p is a viscosity parameter, r p is a relaxation time for 
the polymer conformation tensor, po is the rms extension 
of the polymers in equilibrium, and f(r,t) is a function 
that limits the growth of the trace of R. The model is 
closed by the equation of motion for the conformation 
tensor which reads 



dR. 
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FIG. 1: The mean flow velocity as a function of the distance 
from the wall for the FENE-P (dashed line) vs. the Newtonian 
flow (continuous line). The profiles hardly change near the 
wall, but the amplitude is larger for the FENE-P solution. 
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These equations were simulated on the computer in a 
channel or pipe geometry. The main observations on the 
effect of the polymer on the turbulent flow that we need 
to focus on are the following: (i) For a fixed pressure gra- 
dient at the wall the fluid throughput is increased (see 
Fig. 0). (ii) For a fixed throughput the gradient at the 
wall decreases (i.e. the dissipation decreases), (iii) The 
trace of the conformation tensor R follows qualitatively 
the rms streamwise velocity (see Fig. 2). We are par- 
ticularly interested in point (iii) since in our opinion the 
space dependence of the amount of stretching (and with 
it of the effective viscosity) is crucial, and compare ||, ||] 
for a discussion of this point in the context of the insta- 
bility of laminar flows. Obviously, Eqs. (||)-((|) as they 
stand are not amenable to analytic investigation in the 
turbulent regime. To gain insight we therefore attempt 
to simplify them as much as possible without losing the 
main phenomena (i)-(iii). Consider therefore a model for 
the streamwise velocity which is the Burger's equation 
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FIG. 2: Upper panel: the dependence of the rms velocity 
fluctuations as a function of the distance from the wall. We 
are interested in U Tsns for comparison with our model. Lower 
panel: the trace of the conformation tensor J? as a function of 
the distance from the wall. We stress the qualitative similarity 
to the dependence of U tsnB in the upper panel. 

(it in the streamwise directions with gradients in the y 
(wall-normal) direction), to which the effect of a scalar 
R is added: 

u t + uu y = vuyy + sR y + F , (4) 
1 

Rt + uRy = --R + Ruy. (5) 

where a subscript y stands for a partial derivative with 
respect to y. In the following we shall denote Eqs. (|])- 
with the acronimous uR model. The parameter s is 
related to the polymer concentration, and t is the re- 
laxation time of R. We will consider the model in the 
domain — L < y < L, with boundary conditions chosen 
later. We will denote spatial averages by pointed brack- 
ets, (A) = J A{y)dy. The simplicity of the uR model 
allows us to state the energy budget in simple terms. 
Multiplying (^) by u and taking the spatial average of 
(||) and (H) we obtain: 

\j t {v?)=-v{ul) + s{uR v ) + {Fu) (6) 
±(R) = -±{R)+2(Ru y ) (7) 
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FIG. 3: Comparison of the solution of the uK model (con- 
nected circles) to the solution of Burger's equation (continu- 
ous line) . The dashed line corresponds to the function 7?. The 
parameters are r = 10 4 , uo = 1.17, s = 0.25 and Ro = 2.5. 

The term (uR y ) measures the "energy" given by the ve- 
locity field u to the polymer field R. Multiplying (Q) by 
s/2 and summing (^) with (^]) we obtain: 

J t {\ (u 2 ) + = (u 2 v ) -^(R) + (Fu) (8) 

In the steady state the overall power (Fu) is balanced 
by the overall energy dissipation per unit time D, D — 
v (u 2 y ) + ^p(R). The term \ (u 2 ) + f (R) represents the 
sum of the kinetic energy of the flow plus the potential en- 
ergy of the stretched polymers. We remark that already 
from these elementary consideration it becomes clear that 
for a fixed power input the existence of the term 4= (R) 
necessarily reduces the gradients of u in agreement with 
point (ii) above. To address points (ii) and (iii) further 
we consider next the solution of the model with F = 
and with a fixed velocity u and stretching R at —L and 
L. In other words, we take as boundary conditions 

u(-L) = u , u(L) = -u , R(-L) = R(L) = R . (9) 

In Fig. 3 we compare the solution of the uR model to 
that of the pure Burger's equation (i.e. Eq. dJ) with 
R = F = 0.) To focus our thinking we would like the 
reader to consider the solution in the left half space as a 
model of the streamwise velocity component in the lower 
half channel, with the solution in the right half space 
being simply an anti-symmetric copy. The position of 
the lateral "wall" is modeled by the point where u = 0. 
Thinking this way points (ii) and (iii) are clearly demon- 
strated. We proceed now analytically to demonstrate 
drag reduction (point (ii)) and to understand the profile 
of R (point (iii)). First we consider the stationary solu- 
tion of the pure Burger's equation. Integrate equation 
(Q) in y to find 

iu 2 = vu y + \^u\ (10) 
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where the constant of integration was fixed by noticing 
that for L sufficiently large u y is expected to vanish at 
the boundaries. Multiplying (|T(]) by u y , and integrating 
between —L and L using the boundary conditions, we 
find the viscous dissipation e 

e = v(ul) = |u§ (11) 

Next we consider the solution of the wR model for the 
same boundary conditions (|J) and R > 0. In the station- 
ary state Rt = 0, and by dividing Eq. (||) by Ru we can 
integrate it formally in y and obtain: 

R = a\u\exp(-- [ — ) (12) 

T J U 

where a is a constant of integration. This equation is the 
explanation of point (iii). It says that for small velocity 
u ~ 0, i.e. for y ~ 0, R necessarily goes to 0. In par- 
ticular, approximating u — —my near the point y = 0, 
we obtain R ~ where 6=1 + l/(mr). Thus we 
should expect that at positions with small u where the 
gradient of u is large the generic behavior of R is a cusp 
with R = for mr > 1. To compute the dissipation an- 
alytically we consider the limit r — > oo, i.e. we look for 
a solution at the zero order of the perturbation series in 
1/t. In this limit: 

R=^-\u\. (13) 

1*0 

Returning to Eq. (Q) we integrate it in y to obtain 

vu v = ^u 2 -SR + sRo - ^ul . (14) 

We can now substitute Eq. ( |l3| ) in the domain —L< 
y < where |u| = it, and integrate between — L and 0. 
Multiplying the result by a factor of 2 we find the viscous 
dissipation cr 

2 

£fl = -Uq - 3sRqu (15) 

This result is an analytic demonstration of point (ii). We 
note that our analysis has been performed in the limit 
t — ► oo . For large but finite values of r the qualita- 
tive picture we have drawn is unchanged. Needless to 
say, the above discussion can be reformulated by keep- 
ing constant the energy dissipation while increasing the 
value of u at the boundary, to demonstrate point (i). We 
choose however to demonstrate point (i) next, using a 
forced solution. Point (i) is most clearly demonstrated 
in the wR model using periodic boundary conditions and 
constant forcing. We consider < y < 2tt and choose 
the external forcing F to be: 

F(y) = / sm(4y/3) , for < y < — , 

F(y) = hsin(Ay) , for ^ < y < 2tt . (16) 




FIG. 4: The solution of the uR model with constant forcing 
and periodic boundary conditions (connceted circles). Con- 
tinuous line: the Burger equation without the polymer. 



We examined the solutions of the mR model for the set 
of parameters /o = 0.1, fx = 0.05, v = 0.01, s = 0.01 
and initial conditions R{y) — 2sin(y). The remaining 
parameter is the relaxation time t. It turns out that for 
very small values of r, no effect of the polymer field is 
observed. (We remark that for periodic boundary con- 
ditions the limit r — > corresponds to the case of no 
polymer.) For r — > oo no stationary solutions can be 
obtained. For r smaller than some critical value r c , the 
solution of the itR model shows stable stationary solu- 
tions with drag reduction. The typical situation is pre- 
sented in Fig. [l| showing the numerical solutions for 
t = (0.15) -1 < t c , compared against the solution of 
Burger's equation. The mR model shows a larger ampli- 
tude near the strongest shock due to forcing at x = 3tt/A. 
It is worth noting that the gradient is maintained ex- 
tremely close to the one obtained by the Burger equa- 
tion, demonstrating nicely point (i). Point (iii) is nicely 
demonstrated in Fig. || which presents the solution for R 
together with u 2 for both the uK model and the Burger 
equation. As one can clearly see, the behavior of R is 
similar to what observed in Fig. [l], namely there is a 
qualitative similarity between the space dependence of 
R and it 2 , here with sharp cusp in R near the point of 
maximum gradient of u. On the other hand, the small- 
est shock present in the solution of Burger's equation has 
been completely smoothed out by the uR model. This is 
an indication that when R is not sufficiently suppressed 
where the gradient of u is significant, there can be drag 
enhancement. This important point will be addressed 
again in the concluding remarks. Again, the simplicity 
of the model affords an analytic explanation of why the 
solution near the biggest shock shows a larger velocity 
amplitude compared to the Burger equation. Let yo be 
the position of the maximal velocity near the shock. The 
position yo is unchanged in the two models. We can ex- 
pand u, F and R as power series near yo- Let A be 
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FIG. 5: The solution of the uR model with constant forcing 
and periodic boundary conditions. Dashed line: R(y). Con- 
nected circles: 20u 2 of the solution of uK. Solid line: 20u 2 of 
the solution of Burger's equation. 



denned as A = y s — yo, where y s — 37r/4 is where the 
velocity vanishes (the "wall" position in our model). We 
have: 

u(y) =u - u 2 y 2 , F(y) = F + F x y , R(y) = R + Riy 

07) 

where F ~ / A and Fx ~ — f because of Eq. (16). 
Inserting (0) into the uR model, we obtain: 



u 



vf + sRqt 1 
A/ 



U b (l + P), P 



sRq 



(18) 



where Ub = v / A is the solution obtained without poly- 
mer, (P = 0). The increment of the velocity, which is 
responsible for the drag reduction, is proportional to sRq. 

In summary, we have introduced a simple model of 
the effect of polymeric additives to Newtonian fluids, 
with the aim of understanding in simple mathematical 
terms some of the prominent features associated with 
the phenomenon of drag reduction. Needless to say, the 
model cannot be taken as quantitative; we were con- 
cerned with the qualitative features summarized above 
for convenience as point (i)-(iii). We demonstrated that 
our model reproduces these qualitative features, and pro- 
vided straightforward analytic explanation to all those 
features. It appears that we can draw from the results of 
this model a few important conclusions: (i) Arguments 
concerning the turbulent cascade process do not appear 
essential. These arguments are the hallmark of the theory 
presented in jj] which proposed that the main effect of 
the polymer is to introduce a dissipative cutoff at scales 
larger than the Kolmogorov scale, due to the polymer 
relaxation time matching there the hydrodynamic time 
scale. Looking at Eq. (^) one could think that the largest 
velocity gradients could be estimated as (\u y \) ~ 1/r. 
But since R can go to zero where the gradient is largest, 



such estimates cannot be made. Moreover, just an in- 
crease in the dissipative scale cannot account for drag 
reduction; a homogeneous increase in the effective vis- 
cosity should lead by itself to drag enhancement rather 
than reduction, (ii) Drag reduction is a phenomenon 
that appears on the scales of the system size, involving 
energy containing modes rather than dissipative, small 
scale modes [|l0| . (iii) The main point appears to be the 
space dependence of the stretching of the polymer, here 
modeled by the value of R(y). It is crucial that R is 
small where the velocity gradients are large. It is the 
space dependence of the effective viscosity which should 
be looked at as the source of drag reduction. A similar 
conclusion was arrived to in the context of the study of 
the stability of laminar flows in a channel geometry [|J ||] 
accept that there the space dependence of the effective 
viscosity had been introduced by hand. In the FENE- 
P context as well as in our model (and presumably in 
actual experiments) the space dependence appears self 
consistently. It remains to understand this self consistent 
build up of differential effective viscosity in the context of 
the much more elaborate FENE-P model. In light of the 
present results this appears an extremely worthwhile en- 
deavor that will shed important light on the phenomenon 
of drag reduction. 
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